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INTRODUCTION 


The  author,  as  part  of  her  Joint  Educational  Opportunities  for  Minorities  (JEOM) 
internship  through  the  Department  of  Defense  under  the  cognizance  of  Klaus  Halterman 
of  the  Naval  Air  Warfare  Center  Weapons  Division  (NAWCWD),  China  Lake, 
California,  has  undertaken  learning  computational  physics/electromagnetics  in  order  to 
devise  basic  models  to  analyze  fundamental  problems.  Through  the  use  of  numerical 
methods  and  Matlab,  the  author  modeled  parts  of  problems  that  are  difficult  to  solve 
without  numerical  integration.  This  technical  publication  describes  how  she  used  the 
method  of  moments  as  an  approach  to  solve  integral  equations  in  electrostatic  and 
electrodynamics  problems. 

Line  charges  are  very  fundamental  electrostatic  problems.  By  modeling  the  effect  of 
all  the  charges  in  the  area  surrounding  the  line,  analysts  can  get  an  idea  of  the  electric 
field  around  the  line.  Taking  this  approach  one  step  further,  they  are  able  to  apply  similar 
concepts  to  the  strip  problem.  This  problem  is  set  up  as  a  strip  of  metal  that  has  a 
specific  width,  is  infinitely  long,  and  is  infinitesimally  thin,  a  configuration  that  can  be 
thought  of  as  a  group  of  line  charges  all  lying  next  to  each  other.  By  applying  a  potential 
to  the  strip,  analysts  attempt  to  calculate  the  charge  density  across  the  strip.  In  the 
equations  that  are  used,  the  charge  density  values  are  part  of  the  integral,  thereby  making 
it  necessary  to  use  an  alternate  method  of  numerical  integration.  For  this  task,  the  author 
used  the  method  of  moments.  Herein,  the  strips  are  divided  into  pulses  and  pulse  points, 
and  pulse  matching  is  performed  to  approximate  the  actual  result.  As  such,  the  more 
pulses  that  are  used,  the  more  accurate  the  results  are,  to  a  certain  extent.  In  this  method, 
the  unknown  variable  becomes  a  summation  of  unknown  constants  multiplied  by 
measurable  values  and  can,  therefore,  be  removed  from  the  integral.  By  solving  the 
problem  like  a  set  of  equations  using  matrices,  the  author  was  able  to  find  the  values  of 
the  unknown  constants  and  thereby  detennine  the  charge  density  at  the  pulse  points 
across  the  strip.  The  author  was  able  to  graph  the  charge  density  across  the  strip  and  see 
the  differences  when  different  forcing  functions  were  used,  as  well  as  when  the  number 
of  pulse  points  were  increased. 

The  author  applied  the  same  method  for  alternate  setups  of  the  strip  problems, 
including  coupled  strips  aligned  linearly,  coupled  strip  in  a  cylinder,  and  a  square- 
cylinder  electrodynamics  problem.  The  square-cylinder  problem  was  slightly  different 
because  the  cylinder  was  excited  by  an  incident  plane  wave  nonnal  to  the  cylinder  axis. 
As  a  consequence,  a  transverse  magnetic  scattering  problem  resulted.  It  was  approached 
the  same  way  as  the  other  strip  problems,  but  with  a  different  formula  for  the  forcing 
function.  The  author  was  able  to  look  at  the  current  density  around  the  square  cylinder 
and  at  the  charge  density  in  the  coupled  strip  problems. 
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The  author  examined  a  thin-wire  antenna  problem,  which  she  approached  in  a  very 
similar  way  as  the  strip  problem,  but  used  a  different  set  of  equations.  The  author  was 
able  to  look  at  the  effect  of  different  lengths  of  the  antenna,  different  radii,  and  different 
feed  locations.  She  also  changed  the  forcing  function  from  a  delta  function  excitation  to 
an  incidental  plane  wave,  thereby  generating  a  scattering  problem.  Therein,  she  looked  at 
the  change  in  the  length  of  the  antenna. 


LINE  CHARGES 


In  basic  electrostatics,  a  line  charge  is  a  source  charge  spread  out  along  a  line  of  an 
inconsequentially  small  radius,  which  essentially  acts  as  a  line  of  infinitesimally  small 
point  charges.  The  formula  for  calculating  the  electric  field  due  to  the  line  charges  is 
based  on  the  formula  for  point  charges  (Equation  1). 


E(P) 


4^Sn 


0  <=i  6 


(1) 


This  formula  provides  the  electric  field,  E,  generated  by  the  charges  qi,...,qa  acting 
on  point  P.  where  So  (8.854223  x  10'  Farads/meter)  is  the  permittivity  of  free  space,  r 
is  the  unit  vector  between  the  point  P  and  the  source  charge  c/„  and  r  is  the  distance 
between  the  point  P  and  the  source  charge  q,.  Breaking  the  line  into  infinitely  many 
small  pieces,  each  of  length  dl,  and  calculating  the  charge  at  each  small  piece  result  in  the 
integral  equation  (Equation  2). 


Herein,  r  is  the  unit  vector  between  dl  and  P;  r  is  the  distance  between  dl  and  P;  and  /. 
represents  the  charge  per  unit  length  ( dq/dl ),  which  is  a  constant. 

Figure  1  shows  how  the  results  of  the  equation  change  when  the  distance  varies 
along  the  y-axis.  As  the  data  points  indicate,  an  inverse  relationship  exists  between  the 
distance  for  the  line  and  the  voltage  level.  For  example,  as  the  distance  from  the  line 
decreases,  the  voltage  increases,  and  vice  versa. 

Figure  2  is  a  two-dimensional  plot  of  the  field  strength  around  the  line  charge.  The 
charge  through  the  10-meter  line  is  uniformly  distributed,  and  A  =  5.  Herein,  the  distance 
can  be  varied  in  both  the  x  and  y  directions.  This  figure  shows  the  same  inverse 
relationship  between  the  distance  from  the  line  and  the  voltage  mentioned  earlier. 
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FIGURE  1.  Line  Charge  Problem  Results,  Voltage  vs. 
Distance  From  Line.  The  line  charge  was  centered  at 
x  =  0,  and  the  distances  were  varied  only  in  the 
y  direction. 


FIGURE  2.  Two-dimensional  Plot  of  Field  Strength 
Around  10-meter  Line  Charge. 
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USING  METHOD  OF  MOMENTS  IN  STRIP  PROBLEM 


This  problem  involves  an  infinitesimally  thin,  infinitely  long  strip  of  material  that  is 
of  a  given  width,  centered  along  the  y-axis,  and  has  a  given  potential  applied  across  its 
surface.  The  equation  for  this  problem  is  based  on  the  equation  for  the  potential  of  a  line 
charge  (Equation  3). 


d<f>{p) 


<1l 

2ft  SQ 


In 


P~P 


(3) 


Herein,  d(j)  is  the  electric  potential,  p  is  the  location  at  which  the  potential  is  observed, 
qL  is  the  line  charge  density,  So  again  is  a  constant  (8.854223  x  10'  Farads/meter),  and 
p  is  the  location  of  the  charge. 

The  charge  density  across  the  strip  changes  in  neither  the  y  nor  the  z  direction 
because  of  its  specific  dimensions:  infinitesimally  thin  in  the  z  direction  and  infinitely 
long  in  the  y  direction.  As  a  result,  only  the  x  direction  is  considered  when  calculating 
the  charge  density,  and  the  vector  quantities  p  and  p  may  be  reduced  to  their  scalar 
x-components,  x  and  x  .  The  charge  density  is  found  by  multiplying  the  surface  density, 
qs  (x) ,  by  the  differential  change  in  x  (Equation  4). 

<1l(x)  =  qs{x)dx  (4) 

Placing  this  equation  into  Equation  3  and  integrating  both  sides  results  in  Equation 
5. 


<*>«= f 


■m  qs(x) 


In  t 


x-x 


r  elk 


(5) 


Herein,  ®(x)  is  the  forcing  function  and  8  is  the  width  of  the  strip. 

The  method  of  moments  is  used  to  solve  for  qs  (x)  to  circumvent  the  difficulties  that 
arise  when  solving  for  a  variable  in  an  integrand  in  a  standard  way.  The  first  step  is  to 
divide  the  width  of  the  strip  into  N  divisions  of  SIN  size.  In  this  method,  only  q.  (x)  at 

the  center  of  the  sections  is  measured  and  that  measurement  is  utilized  for  the  whole 
section  (see  Figure  3).  This  method  involves  using  pulse  expansion,  wherein  the  pulse 
function  is  integrated  into  the  integral  through  a  summation. 
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FIGURE  3.  Pulse  Expansion. 

The  pulse  expansion  function  is  defined  as  Equation  6. 

q{x)  =  qs(x)*YjI»P"(x)  (6) 

n= 1 

Herein,  q(x)  denotes  the  true  surface  charge  density,  q.  (x)  designates  the  unknown 

function,  N  is  the  number  of  sections  into  which  the  function  is  split,  Ii  through  In 
represent  unknown  constants,  and  pi(x)  through  Pn(x)  are  assumed  known  functions.  The 
summation  is  approximately  equal  to  qs  (x) .  Each  pn(x)  is  a  pulse  function,  as  defined  in 
Equation  7. 


1, 

0, 


x„ 


S 

2  N 


<x<xn  + 


S 

2  N 


otherwise 


Herein,  xn  denotes  the  nth  pulse  point. 


(7) 


When  this  pulse  function  is  incorporated  into  the  integral  equation  (Equation  5)  for 
anything  beyond  the  limits  of  x„  -  ( S/2 )  and  xn  +  (<572)  [xn  -  (SUN)  and  x„  +  (SUN)),  the 
function  disappears;  so,  the  integral  is  represented  by  Equation  8. 


(8) 
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Next  is  the  testing  stage,  in  which  the  summation  is  expanded  with  different  n’s, 
depending  upon  the  number  of  sections  into  which  the  function  is  split,  a  value 
previously  referred  to  as  N.  By  setting  n  =  1,2,3  analysts  can  create  an  Ax  A  matrix 
a  in  which  each  row  (m)  represents  the  potential  on  a  section  and  each  column  (n) 
represents  how  that  section  affects  the  pulse  point  at  xm.  In  the  matrix  element  amn,  if  H 
represents  the  width  of  the  pulse  and  w  represents  the  width  of  the  strip,  then  a„  =  (-w/2) 
+  (n  -  \)H  is  the  lower  limit  of  the  pulse,  b„  =  (-w/2)  +  nH  is  the  upper  limit  of  the  pulse, 
and  xm  =  (-w/2)  +  (n  -  Vi)H,  which  puts  the  points  in  the  center  of  the  pulses.  The 
following  functions  (Equations  9  and  10)  are  used  to  calculate  the  unknown  values. 


«11  • 

a\N 

"Al 

_aN\ 

aNN  _ 

1 - 

1 _ 

1 

•  ^ 

_ 1 

a 


mn 


1 

2;rs0 


L  =  potential  at  xm  =  0(x„, ) 


(9) 


(10) 


By  multiplying  the  vector  /  by  the  inverse  of  the  matrix  a,  one  can  solve  for  the 
unknown  I  ,/s.  Problems  arise  when  calculating  the  integral  of  the  diagonal  terms  (an) 
due  to  the  singularity  at  xm,  where  the  denominator  in  the  integral  turns  to  0.  Gaussian 
quadrature  is  used  to  avoid  this  singularity  by  splitting  the  pulse  in  half  and  summing  the 
integrals  from  both  sides.  This  type  of  quadrature  does  not  evaluate  the  integrand  at  the 
end  points,  thereby  avoiding  this  singularity. 

The  graphs  that  follow  (Figures  4  through  8)  represent  the  results  of  this  problem 
solved  via  the  numerical  method  described.  Figures  4,  5,  and  6  show  the  results  of 
solving  the  problem  with  100  pulses.  Figure  4  shows  the  data  for  this  problem  when  the 
forcing  function  is  &(x)  =  1.  In  contrast,  Figure  5  shows  the  results  when  the  forcing 
function  is  <P(x)  =  x,  an  odd  function.  By  comparing  the  two,  analysts  see  that  the 
function  is  odd  when  the  forcing  function  is  odd  and  is  even  when  the  forcing  function  is 
even.  This  outcome  is  confirmed  by  Figure  6,  which  shows  the  results  for  a  forcing 
function  of  &(x)  =  x  . 

Figures  7  and  8  afford  a  comparison  of  the  results  for  similar  forcing  functions  with 
different  numbers  of  pulses.  The  curve  shown  in  Figure  7  represents  the  data  for 
15  pulses  and  that  depicted  in  Figure  8  provides  the  data  for  50  pulses.  These  figures 
show  that  the  function  becomes  much  smoother  as  more  pulses  are  used;  although,  the 
curve  does  not  become  significantly  smoother  as  the  number  of  pulses  increases 
beyond  50. 
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FIGURE  4.  Strip  Problem  Results,  Charge  Density 
Across  2-meter-wide  Strip,  With  Even  Forcing  Function 
of  =  1 . 


FIGURE  5.  Strip  Problem  Results,  Charge  Density 
Across  2-meter-wide  Strip,  With  Odd  Forcing  Function 
of  Ofx)  =  x. 


NAWCWD  TP  8673 


FIGURE  6.  Strip  Problem  Results,  Charge  Density 
Across  2-meter-wide  Strip,  With  Even  Forcing 
Function  of  <&(x)  =  x  . 


FIGURE  7.  Strip  Problem  Results,  Charge  Density 
Across  2-meter-wide  Strip.  Herein,  the  number  of 
pulses  (AO  is  15  and  the  forcing  function  is  <b(x)  =  1. 
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Strip  Problem:  Solving  for  Charge /Unit  Length  given  width=2,  §(<e)  =  i 


FIGURE  8.  Strip  Problem  Results,  Charge  Density 
Across  2-meter-wide  Strip.  Herein,  the  number  of 
pulses  (AO  is  50  and  the  forcing  function  is  <5>(x)  =  1. 

Another  variation  on  the  strip  problem  involves  incorporating  two  strips  (coupled 
strips)  into  the  equation  and  adjusting  the  calculations  to  reflect  the  presence  of  another 
strip.  This  addition  is  first  examined  as  the  two  strips  are  aligned  along  the  same  axis  and 
are  separated  only  by  a  small  amount  of  space.  In  these  calculations,  the  matrix  becomes 
a  partitioned  matrix  and  is  only  symmetric  within  each  portion.  The  matrix  is  set  up  as 
follows. 

effects  of  A  on  A  effects  of  B  on  A 
effects  of  A  on  B  effects  of  B  on  B 

Figure  9  shows  the  outcome  of  applying  this  method  to  coupled  strips  that  are 
aligned  along  the  same  axis.  The  figure  indicates  that,  though  the  strips  are  not 
connected,  the  existence  of  another  set  of  charges  influences  the  charge  distribution. 

In  another  configuration  used,  the  strips  were  arranged  in  a  cylinder  and  infinitely 
small  gaps  where  left  where  the  strips  would  meet.  The  data  in  Figure  10  reflect  the 
outcome  of  this  arrangement,  and  the  figure  shows  that  the  singularity  that  occurs  at  the 
edge  of  the  strip  is  still  present,  though  the  direction  of  the  singularity  is  opposite  from 
strip  to  strip.  This  trend  is  also  clear  in  Figure  1 1;  the  data  in  that  figure,  which  represent 
the  charge  density  measured  around  the  cylinder,  is  plotted  in  two  dimensions — radians 
and  charge  density. 
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FIGURE  9.  Coupled  Strip  Problem  Results,  Charge 
Density  Along  Coupled  1 -meter- wide  Strips  Aligned  on 
X-axis  and  Separated  by  0.1 -meter- wide  Space. 


Density  Around  Coupled  Strips  Oriented  in  Cylindrical 
Shape  With  1 -meter  Radius. 
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FIGURE  11.  Two-dimensional  Plot  of  Charge 
Density  Around  Coupled  Strips  Oriented  in 
Cylindrical  Shape  With  1 -meter  Radius. 


Yet  another  variation  on  this  problem  is  the  square  cylinder  scattering  problem  in 
which  one  strip  is  shaped  into  a  square-cylinder  and  excited  by  an  incidental  plane  wave. 
Through  orientation  of  the  cylinder  in  such  a  way  that  the  plane  wave  is  propagating  in  a 
direction  normal  to  the  axis  of  the  cylinder,  the  induced  surface  current  and  scattered 
electric  field  have  components  only  in  the  direction  of  the  axis  of  the  cylinder.  This 
situation  is  called  transverse  magnetic  illumination  and  is  solved  with  an  electric  field 
integral  equation  (EFIE).  Because  this  is  a  two-dimensional  problem,  the  EFIE  is 
represented  by  Equation  1 1 . 

K(p)  =  ^\j,(p')nfik\p-p'\)dl',  psC  (11) 

Herein,  e[  =  designates  the  z-component  of  the  incidental  electric  field,  p  is  in  radians 
(\p-  p'\  ( p '  denotes  the  source  point)  is  the  distance  between  the  points,  k  is  the 

wavenumber  (2tc/A,  in  which  /.  is  the  wavelength),  rj  is  the  intrinsic  impedance  (377  ohms 
for  air),  C  designates  an  arbitrary  closed  contour,  Jz  is  the  induced  current  density 
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(amperes/meter),  and  H „2)  is  a  zero-order  Hankel  function  of  the  second  kind.  The 
distance  portion  can  be  rewritten  in  Cartesian  coordinates  as  Equation  12. 


|  P~P' 


(12) 


Herein,  xn,  yn,  xm,  and  ym  denote  Cartesian  coordinates  that  represent  the  intervals  along 
the  circumference,  x  designates  the  unit  vector  in  the  x  direction,  y  designates  the  unit 

vector  in  the  y  direction,  /'  denotes  an  arc  length  variable,  and  ln  represents  a  unit  vector 
directed  along  the  nth  segment. 

The  diagonal  tenns  are  not  as  easy  to  handle  as  in  the  previous  method  of  moments 
cases;  so,  instead,  the  fonnula  inside  the  integral  is  rewritten  as  Equation  13  to  avoid  any 
problems. 


H, 


(2)1 


l/'l)= 

ln 

'kr' 

+r 

+i-j- 

ln 

'kV' 

+  r 

n 

l  2  J 

n 

l  2  J 

(13) 


2 

Herein,  /  designates  length;  j  =  -1;  and  y  denotes  Euler’s  constant,  approximately  0.5772. 

The  final  factor  to  consider  is  the  “forcing  function”  or,  here,  an  incidental  plane 
wave,  which  is  calculated  using  Equation  14. 

El  (p)  =  /•;'e/AUc‘“''-v'i:''y  1  (14) 


Herein,  9‘  is  the  angle  of  the  incidental  plane  wave  as  measured  from  the  x-axis. 

Solving  for  the  current  by  using  the  method  of  moments,  analysts  can  determine  the 
distribution  of  current  around  the  cylinder  (see  Figure  12).  This  figure  provides  the 
results  for  only  half  way  around  the  cylinder,  starting  at  the  middle  of  the  side  facing  the 
incidental  wave.  The  data  indicate  that  the  current  is  concentrated  on  the  side  that  comes 
into  contact  with  the  incident  wave  first.  As  the  reader  can  see,  the  singularities  occur  at 
the  comers  (the  corners  are  at  pulses  10  and  30).  Also  apparent  is  the  lack  of  current  on 
the  far  side. 
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Scattering  Problem:  Square  Cylinder,  a=l,  A=i,  N=80 


FIGURE  12.  Square-cylinder  Scattering  Problem  Results, 
Normalized  Current  Density  Half  Way  Around  Square 
Cylinder  Excited  by  Plane  Wave  Oriented  at  Angle  of 
71  Radians  From  X-axis. 


THIN- WIRE  ANTENNA  PROBLEM 


A  wire  is  deemed  thin  when,  due  to  the  thinness  of  the  wire,  one  can  assume  the 
surface  current  ( J  )  is  the  sum  of  the  interior  and  exterior  surfaces  of  the  cylinder.  This 
assumption  is  allowable  when  the  thickness  of  the  wire  is  small  in  comparison  with  the 
wavelength  of  the  incidental  wave  (/.),  typically  around  0.01k.  Because  the  wire  is  so 

thin,  E‘  is  invariant  around  the  cylinder,  thereby  leaving  the  only  component  of  E‘ 
tangential  to  the  surface  to  be  the  component  in  the  z  direction.  In  this  problem,  an 
inherent  assumption  is  that  the  antenna  is  a  perfect  electrical  conductor,  thus  allowing  the 
boundary  condition  represented  in  Equation  15. 

El  +  £'  =  0  (15) 

Herein,  Esz  is  the  z-component  of  the  scattered  electromagnetic  field,  and  E\  is  the 
z-component  of  the  incidental  electromagnetic  field.  The  antenna  is  also  said  to  have  a 
radius  of  a  and  a  height  of  2 h.  Through  derivation,  one  can  generate  Equations  16 
and  17. 
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k2\  k2A,  +- 


Az=^f  Jz(z)f 
z  4 ni--=-h  zV  ’ie 


-a  d6  dz , 


Herein,  co  =  frequency;  R  =  r  -  r '  ;  k  is  the  wavenumber,  defined  as  k  =  co^f/us  ,  where  // 

=  the  permeability  of  the  medium  and  s  is  the  permittivity  of  the  medium;  Az  =  the  vector 
potential;  and  / 2  =  -1 . 

To  solve  Az,  Equations  18,  19,  and  20  are  used.  Herein,  z  is  a  variable  used  to 
describe  a  coordinate. 


/  .  \  1  f  n  £  j  R 

K[z-z  )  =  —  f  - - dO 

v  ’  2n]-*  R 


Jz(z)  =  2  7ia(l(z)) 


lf/(z)4z-z)rfz 


4;r£  Jl  <3z 


If  /(z )  =  J  /(z )  k(z  -  z)dz  ,  then  Equation  2 1  applies. 


i  „  V/52  3 

E,^--\Mh+ky^ 


The  solution  to  this  differential  equation  is  in  the  form  of  Equation  22. 

f(x)  =  fh(x)  +  fP(x) 

Herein,  f,  is  the  homogeneous  solution  and fp  is  the  particular  solution. 


Solving  as  a  differential  equation  results  in  Hallen’s  equation  (Equation  23). 
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EUZ)  =  J 


V 


Ank 


f  d2 
~d? 


+  k 1  |  /(z)if  (z-z)Jz 


=  C cos kz  +  B sin kz  +  v(z),z  e  ( -h,h ) 


v(z)  =  — J  fi'(z)sinA(  |z-(z|  ))dz  ,  with  /(±/j)  =  0  (23) 

k  ^ 

This  equation  can  be  solved  by  using  the  method  of  moments.  For  an  antenna  problem, 
E[  is  replaced  by  VS  ( z  -  zg  j  (the  Dirac  delta  function),  V  denotes  the  delta  function 

“strength,”  and  the  equation  is  solved  by  using  pulse  expansion  and  testing  for  a  current 
fed  at  zg.  Based  on  the  assumption  that  the  forcing  function  is  fed  at  z0,  the  delta  function 
is  evaluated  at  zq.  The  first  step  is  to  expand  with  pulses  of  uniform  width  A  and 
approximate  I(z)  via  Equation  24. 


h4=27„n»6) 


I,'  s  are  unknown  complex  constants  and  Equation  25  applies. 


nH(z)  = 


.  A  A 

!,  tf  zn--<z<zn+- 

1,  otherwise 


(24) 


(25) 


Half  of  a  pulse  is  placed  at  either  end  of  the  antenna  and  given  a  height  of  zero  to 
account  for  the  constraint  of  0  current  at  +h  and  -h.  Then,  enforcing  the  resulting 
equation  at  the  matching  points  (zm)  results  in  Equation  26. 


(zm )  =  j[  Tt\  f /,  X  7»n»  (iK  (z,n  ~z)dz-C  (cos  kzm  )-B(  sin  kzm ) 


.  4  7ik 

Herein,  m  =  0, 1 ,2, . . . N  +  1 . 


(26) 


This  equation  can  be  rewritten  when  considering  the  delta  function  as  Equations  27 
and  28  to  derive  Equations  29,  30,  and  3 1 . 


F 


=  vm(z)  =  lr^(z)sm,(k 


(27) 


jy 

V  /  Z  -  C  cos  Az  -  5  cos  kz  =  V  ,  m  =  0, 1, 2, . . .  N  + 1 

72  22272  722  722  722  ‘  ‘  ’  ’ 


72=1 


(28) 
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zm  =  -h  +  mA,  Zmn  =  j  1  "J  K  (zm  -  z)  dz, 


(29) 


1  -y’M 

K(z-z)  =  —  [  - dd,  R  =  \r  -r'\  =  .  (z-z)2  +4 a2  sin2  — , 

V  ’  2n]-*  R  1  1  VV  ’  2 


(30) 


A  = 


2  h 

N  +  l 


(31) 


There  are  now  N  +  2  equations  (including  the  pulses  at  either  end)  and  A+2  unknowns 
(including  B  and  C),  which  can  be  solved  using  linear  algebra.  By  setting  up  the 
following  matrices,  one  can  solve  for  the  unknown  F s,  as  well  as  C  and  B. 

cos  (ykzi  0 )  &  z^  0 1  & . . .  &  z^  ON  &  -  sin  ( kzi  0)  @-  cos  [kz^  l)  &  z±  1 1  & . . .  &  Zl  1 N  & 
-sin(kzil)@:&:&'-.&:&:@-  cos(kziA)  &  z±  (A^,l)  .  &  Z^NN  &  -sin  (kz;iV)@-  cos(kz;  («  + 1)^) 


The  rows  (m)  each  represent  a  point  on  the  wire,  while  the  columns  (n)  each 
represent  the  corresponding  pulse  on  the  wire  and  the  way  in  which  the  pulse  affects  the 
other  points. 

The  graphs  provided  in  Figures  13  and  14  show  the  different  variations  of  this 
problem.  Plotted  in  Figure  13  is  the  magnitude  of  the  complex  vector  when  2=1  meter, 
N=  53,  a  =  0.007022  meter,  and  the  length  of  the  rod  is  A.  Figure  14  shows  the  change 
in  the  current  when  the  size  of  the  radius  ( a )  is  decreased. 

Figure  15  shows  the  distribution  of  current  changes  as  the  location  of  feed  input, 
which  was  previously  held  at  z  =  0,  is  moved.  These  data  indicate  how  the  current 
changes  as  it  is  moved  repeatedly  three  pulses  away.  Here,  the  asterisks  (*)  represent  the 
location  of  input,  the  dotted  lines  designate  the  imaginary  parts,  and  the  solid  lines 
signify  the  real  parts.  Figure  16  shows  the  plot  of  the  magnitude  of  the  vector. 
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FIGURE  13.  Antenna  Problem  Results,  Current 
Along  Thin- wire  Antenna  (Length  =  A/2). 


FIGURE  14.  Antenna  Problem  Data  Showing 
Differences  Resulting  From  Changes  in  Size  of 
Radius  (a),  Current  Along  Thin-wire  Antenna 
(Length  =  A). 
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FIGURE  15.  Antenna  Problem  Results  Showing  Effect  of 
Varying  Input  Location  of  Current  Along  Thin-wire 
Antenna. 


FIGURE  16.  Antenna  Problem  Results,  Magnitude  of 
Current  Vector.  Graph  shows  effect  of  varying  current  input 
location  along  thin-wire  antenna. 
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The  graphs  in  Figure  17  provide  a  comparison  of  the  effect  that  different  lengths  of 
wire  have  on  the  resultant  data.  These  graphs,  plots  of  the  magnitude  of  the  complex 
vector,  show  the  sine  curves  for  the  input  force.  As  the  reader  can  see,  the  current  of  this 
wire  is  dependent  upon  its  length. 


FIGURE  17.  Antenna  Problem  Results  for  Different  Lengths  of  Thin- wire  Antenna. 
Graphs  provide  a  comparison  to  demonstrate  the  effect  on  the  current 
by  varying  the  thin- wire  antenna  length  (0.5,  1,  1.5,  and  2  meters). 


Changing  the  force  acting  upon  the  antenna  to  an  incidental  wave  results  in  the 
problem  becoming  a  scattering  problem.  Figure  18  shows  the  distribution  when  the  wire 
length  is  A./4,  and  Figure  19  shows  the  scattering  distribution  when  the  wire  is  k!2. 
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FIGURE  18.  Scattering  Problem  Results,  Current  Along 
Thin-wire  Antenna  (Length  =  A./4)  Excited  by  Incidental 
Plane  Wave. 


FIGURE  19.  Scattering  Problem  Results,  Current  Along 
Thin-wire  Antenna  (Length  =  X/2)  Excited  by  Incidental 
Plane  Wave. 
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CONCLUSION 


This  document  describes  the  author’s  experience  as  an  intern  at  the  NAWCWD  in 
learning  the  basics  of  electrostatics  and  electrodynamics,  as  well  as  how  to  approach 
computational  models  in  physics.  Herein,  she  discusses  the  models  and  methods  that  she 
adopted  to  solve  various  fundamental  problems  applicable  to  these  fields.  Moreover,  the 
author  shows  how  to  use  the  method  of  moments  to  calculate  the  charge  density  across  a 
strip  and  the  current  along  a  thin- wire  antenna. 


NOMENCLATURE 


12 

So  permittivity  of  free  space  (8.854223  x  10'  "  Farads/meter) 
s  permittivity  of  a  medium 

p  permeability  of  a  medium 

A,  charge  per  unit  length 

wavelength 

(f>  electric  potential 

0(x)  forcing  function 

S  width  of  strip  in  strip  problem 
p  distance  between  the  points 

p-  p'\  P  in  radians 

i/  intrinsic  impedance 

y  Euler’s  constant,  approximately  0.5772 

Q'  angle  of  incidental  plane  wave  as  measured  from  x-axis 
ciu  fth  diagonal  term  of  matrix  a 

amn  (m,n)’th  element  of  matrix  a 

a„  lower  limit  of  pulse  n 

b„  upper  limit  of  pulse  n 

C  an  arbitrary  closed  contour 

E  electric  field 

H  width  of  pulse 

H\^  zero-order  Hankel  function  of  the  second  kind 
In  unknown  constants 

j  the  square  root  of  -1,  i.e.  the  imaginary  unit 

j  surface  current 

Jz  induced  current  density  (amperes/meter) 

k  wave  number  (2n/X,  in  which  /.  is  the  wavelength) 
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/ 

K 

In 

N 

P 

P 

P' 

P 

P 

Pn(x) 

q 

p 

q»  •••»  qn 
qL 
q(x) 

qs  (*) 


NAWCWD 

V 

w 

X 


X 

X 
X  iii 
X„ 

y 

y 

y>n 

yn 

Z 

z 


length 

an  arc  length  variable 

a  unit  vector  directed  along  the  nth  segment 
number  of  pulses  used  in  a  method  of  moments  solution 
a  point  at  which  electric  field  strength  is  observed 
an  arbitrary  point  along  C 
source  point 

location  at  which  potential  is  observed 
location  of  a  charge 

pulse  functions 
electric  charge 
source  charge 
charges 

line  charge  density 
true  surface  charge  density 
potential  across  surface 

distance  between  the  point  P  and  source  charge  q,  (Equation  1) 
distance  between  dl  and  P  (Equation  2) 
unit  vector  between  point  P  and  source  charge  q,  (Equation  1) 
unit  vector  between  dl  and  P  (Equation  2) 

Naval  Air  Warfare  Center  Weapons  Division 
Dirac  delta  function  “strength” 
width  of  strip 

horizontal  direction  in  two-dimensional  space 

one  of  the  three  fundamental  directions  in  three-dimensional  space 

x-component  of  p 

x-component  of  p 
unit  vector  in  the  x  direction 
x-coordinate  of  an  arbitrary  point  m  along  C 
x-coordinate  of  an  arbitrary  point  n  along  C 
vertical  direction  in  two-dimensional  space 

one  of  the  three  fundamental  directions  in  three-dimensional  space 

unit  vector  in  the  y  direction 

y-coordinate  of  an  arbitrary  point  m  along  C 

y-coordinate  of  an  arbitrary  point  n  along  C 

one  of  the  three  fundamental  directions  in  three-dimensional  space 

variable  used  to  describe  a  coordinate 
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